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One-dimensional Bose gas with contact interaction in optical lattices at zero temperature is in¬ 
vestigated by means of the exact diffusion Monte Carlo algorithm. The results obtained from the 
fundamental continuous model are compared with those obtained from the lattice (discrete) Bose- 
Hubbard model, using exact diagonalization, and from the quantum sine-Gordon model. We map 
out the complete phase diagram of the continuous model and determine the regions of applicability 
of the Bose-Hubbard model. Various physical quantities characterizing the systems are calculated 
and it is demonstrated that the sine-Gordon model used for shallow lattices is inaccurate. 


The Bose-Hubbard model (BHM) was introduced in 
1963 [H 12]. While the original motivation was to de¬ 
scribe a crystalline solid, for which the model failed, the 
BHM became one of the fundamental quantum many- 
body problems. It has found clear-cut realization with 
ultracold atoms in deep optical lattices. This lead to the 
seminal observation [3] of the superfluid-Mott-insulator 
quantum phase transition [3] following the proposal of 
Ref. [Sj. In many aspects the experiments surpass the 
theory as shallow optical lattices can be easily realized, 
while no exact quantum many-body description of such 
systems is known up to date. Even the case of deep op¬ 
tical lattices is controversial, as the scattered discussions 
demonstrate |6UTT] indicating the necessity to go beyond 
the standard BHM (for a review see [12] )• Neverthe¬ 
less, the BHM is commonly used for lattice systems in 
different dimensions and it frequently works very well. 
Still, there arise natural and important questions that 
motivated the present work: When can it be used with 
confidence? What is the regime of validity of the BHM? 

The discrete BHM is derived from a continuous space 
model that, due to its complexity, has only been ad¬ 
dressed recently [T^HTS] . In this Letter, we use the ex¬ 
act diffusion quantum Monte Carlo method naimiii- 
m and investigate one-dimensional Bose gas in optical 
lattices using a continuous Hamiltonian in real space. 
We compare the results with those obtained from the 
BHM and determine its regions of validity. Furthermore, 
a whole new generation of clean experiments on one¬ 
dimensional Bose gases loaded in optical lattices [221 [ 25 ] 
have appeared, while the comparison of theory with ex¬ 
periment is not perfect [211 [25]. We also analyze the 
sine-Gordon (SG) model, commonly used for shallow lat¬ 
tices 0121117 ], and show that it cannot be straightfor¬ 
wardly used for predicting for the position of the phase 
transition and the value of the gap [ 21 ]. We calculate 
the static structure factor, the one-body density matrix, 
the energy gap, the Luttinger parameter and its depen¬ 
dence on the interaction and lattice strengths. Finally, 
we compare our results with the experiments of Ref. [24j 
- surprisingly our theory, which is in principle superior 


to all approximate ones, does not always provide a better 
description. 

The first quantization Hamiltonian of N bosons of 
mass m interacting by a contact potential of strength 
g = —2lT?/(maiB), with oid being the one-dimensional 
s-wave scattering length, has the form 
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The external potential Vi,{x) = Vb cos^(7ra:/ao) repre¬ 
sents an optical lattice of strength Vg with the lattice 
constant oq. A characteristic energy associated with the 
lattice is the recoil energy Free = TT^f?/{2mal). We con¬ 
sider a system of finite size Lao, where L is an integer, 
and impose periodic boundary conditions. 

The ground-state properties of Hamiltonian (1^ are 
studied using the Diffusion Monte Carlo (DMC) algo¬ 
rithm [28j that solves the Schrodinger equation in imagi¬ 
nary time. Statistical variance is significantly diminished 
by using the importance sampling. The DMC method 
gives exact estimation of any observable commuting with 
the Hamiltonian, and delivers bias-free predictions for 
other observables by pure estimator techniques [28] . 

In deep optical lattices, model Eq. 0 reduces to the 
BHM. In its standard and simplest form, the second 
quantization Hamiltonian is given by 
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The hopping and interaction constants J and U are de¬ 
termined as m 

pLao 
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where We{x) is the Wannier function for the lowest Bloch 
band (maximally) localized near the minimum x = of 
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the periodic potential Vi^{x). The results for the BHM 
are obtained by exact diagonalization. 

Figure presents the complete phase diagram of the 
continuous model compared with various theories. The 
transition line separating the superfluid and the Mott in¬ 
sulator phases is obtained from the Luttinger parameter 
K = v-p!c, where the Fermi velocity t;F = hnNKLaom) 
is entirely fixed by the system setup, while the speed of 
sound c depends in a non-trivial way on the strength of 
the interaction and the lattice height [HIHZ]- For unit fill¬ 
ing, the phase transition takes place at the critical value 
iF = 2 as it follows from effective renormalization group 
theory m- The phase transition of the continuous model 
is shown by black circles. For Vb = 0 the critical value 
IoidI/oo = 0.56 coincides with that of the Lieb-Liniger 
model. It is interesting to compare with the SG and 
BHM, which are expected to be valid for shallow lattices 
and high lattices with weak interactions, respectively. As 
there is no way to establish the exact regions of applica¬ 
bility of each theory internally, we deduce them by direct 
comparison with the DMC data. Within the BHM, the 
transition is governed by a single parameter, J/U. Us¬ 
ing the K = 2 criteria for = 12 particles one obtains 
a critical point at {J/U)c = 0.28. The relation to the 
two parameters of the continuous model (the lattice in¬ 
tensity and the interaction strength) is obtained from 
Eqs. (^, resulting in the solid green line in Fig. We 
find that for Vb/Uj-ec ^ 3 the BHM and the continuous 
model predict the same transition curve. There is a cer¬ 
tain discrepancy between the two models at lower ratios, 
as for instance one gets IoidI/uo = 0.402 for = 12 and 
IoidI/oo = (4/3)(J/C/)c = 0.373 in the thermodynamic 
limit for Vb = 0 in the BHM. However, deviations are 
not as dramatic as for other quantities, for instance the 
one-body density matrix which is discussed later. The 
SG model should be valid for shallow lattices, but it is 
not clear, within the model, up to which maximum value 
of Vb it works. It is kind of a surprise that the SG model 
coincides with the DMG result only at Vb = 0, deviating 
from it for any finite value of Vq. There is an overall 
good agreement with the experimental position of the 
phase transition [24]. In the region of shallow lattices, 
the amplitude modulation measurements are compatible 
with the DMG results, while the transport measurement 
at the weakest lattice agrees better with the SG model. 

Figure [^reports the Luttinger parameter K of the con¬ 
tinuous model as a function of |aiD|/ao for a number of 
characteristic values of Vq. The figure also shows the 
BHM prediction and the Vq = 0 Lieb-Liniger limit. Its 
knowledge is essential in order to use the Luttinger liquid 
(LL) theory, which provides a description of long-range 
and small-momentum correlation functions. It is impor¬ 
tant to realize that the effective LL theory uses K as an 
input, while a full quantum many-body problem needs to 
be solved in order to obtain the dependence of K on the 
system parameters. The condition K = 2 provides the 
critical value of |aiD|/ao corresponding to the superfluid- 
insulator transition. For Vb = 0 the line starts exactly 



FIG. 1: (Color online) The zero-temperature phase diagram 
of the continnous model compared with various approximate 
theories, as a function of the s-wave scattering length aio /fflo 
and the optical lattice intensity Vo/Erec- The position of 
the Mott-insulator-superfluid phase transition is defined by 
K = 2. Black circles, DMC results; green solid line, exact 
diagonalization of the BHM; dashed red line, sine-Gordon 
predictions; squares with error bars, (pink - amplitude mod¬ 
ulation, blue - transport measurements) experimental results 
of Ref. |24| . DMC and exact diagonalization results are ob¬ 
tained for N = 12. 

at AT = 1 (for Tonks-Girardeau Cid = 0 gas) and in¬ 
creases with the scattering length. In that case the DMG 
results are compatible with the Bethe ansatz solution of 
the Lieb-Liniger model in the thermodynamic limit, while 
the deviations at weak interactions can be attributed to 
finite-size corrections. In the Mott insulator regime, the 
sound is absent resulting in vertical lines for AT < 2 in the 
thermodynamic limit [29] . Within the BHM, K depends 
on the single parameter J/U, generating a series of curves 
scaled by the value of aio- For shallow lattices BHM pre¬ 
dictions lie above the Lieb-Liniger Vb = 0 curve, which 
by itself can serve as a test of validity of the BHM. A 
more precise boundary of applicability is obtained when 
compared with the DMC results. We find good agree¬ 
ment for large |aiD|/ao and Vb. For the range reported 
in Fig. [^agreement for AT is achieved for Vb/Arec ^ 4. 

The static structure factor [5DH55] S{k) is defined as 
‘S'(fc) = ^{^p{k)^p{—k)) where Ap{k) is the Fourier 
transform of the density-fluctuation operator (following 
the notations from Ref. m)- From the solution of the 
BHM it can be obtained as m 

Sik) « SBKik) = 1 + Gl{k) [So{k) - 1] , (4) 

where Go(fc) = dx \ Wi{x)f' exp [—ik [x — x^)] and 

So{k) = {{ni^ni^)-{ni,){ni 2 ))exp [ikao (h - U)](5) 

is a discrete analogue of S{k). The typical behavior of 
S{k) is shown in Fig. In the Lieb-Liniger gas, corre- 
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FIG. 2: (Color online) Parameter K — 2S'(fcmin)fcF/femin cal¬ 
culated at the first point, fcmin = 27 r/(Lao), of the static 
structure factor for = 12 particles as a function of in¬ 
teraction strength |aiD|/a for lattices of different heights 
Vo = 0,1, 2, 3,4, 5i?rec (top to bottom curves). Solid symbols 
connected with solid lines, DMC results; open symbols, BHM; 
dashed line, Lieb-Liniger thermodynamic result for Vb = 0. 
The K = 2 short-dashed line separates the superfluid regime 
(above) from the insulating one (below). Notice that K can 
be identified with the Luttinger parameter in the region with 
K >2. 


spending to Vb = 0, and large oid the S{k) is a feature¬ 
less monotonous function typical for weakly interacting 
Bose gas [Ml |3S]. The limit of Vq = 0 and Oid = 0 
corresponds to Tonks-Girardeau gas [3S] which can be 
mapped to an ideal Fermi gas showing a kink at fc = 2kp, 
with the Fermi momentum kp = ttN/{L ao). For a fi¬ 
nite optical lattice, more features appear at momentum 
k = 2k-p, which corresponds to the border of the first 
Brillouin zone. 

At low momenta, the static structure factor is well ap¬ 
proximated by the Feynman relation: 

S{k) = 2^ , e{k) = + [hekf , (6) 

where A > 0 is the energy gap. Note that in units oik/k-p 
there is a critical value of the slope in S(k) corresponding 
to AT = 2, see Fig. If S(k) lies above the critical line 
for small momentum, the gas is superfluid and S{k) is 
linear for /c —>■ 0, otherwise the system is insulating and 
S{k) is quadratic for k ^ 0. 

The coherence properties differ significantly in the in¬ 
sulating and the superfluid phases. Commonly, the su¬ 
perfluidity is associated with the presence of a Bose- 
Einstein condensate, where the condensate wave func¬ 
tion is the order parameter of the superfluid phase, man¬ 
ifested by the off-diagonal long-range order (ODLRO) 
in the one-body density matrix (OBDM) pi(r). A fi¬ 
nite valne of the condensate fraction, pi(r —)■ oo)/p ^ 0, 
was used to localize the superfluid-Mott-insulator phase 



FIG. 3: (Golor online) Static structure factor S{k). Solid 
symbols, continuous model; solid lines, SBii{k) defined by Q; 
dashed line, So(fc) defined by dash-dotted line, S{k) = 
|fc|/fcF corresponds to a linear slope with the critical value of 
the Luttinger parameter K — 2 and separates superfluid (SF) 
and Mott-insulator (MI) phases. The following parameters 
are used: short-dashed line, Tonks-Girardeau gas (i.e. Vb = 0 
and fliD = 0); continuous and Bose-Hubbard models (from 
top to bottom) Vo/Erec = 0;4;8 and |aiD| = 4ao. 


transition in three dimensions [TB]. However, in one di¬ 
mension quantum fluctuations destroy ODLRO, even at 
zero temperature [37], and pi(r) always decays to zero. 
The LL theory predicts a slow power-law decay in the 
superfluid phase [BB], in contrast with the fast exponen¬ 
tial decay in the insulating phase. Figure [^ shows the 
OBDM pi{r) calculated in two different phases at the 
largest possible length r = L/2 in a box of size L. The 
solid lines show fits to the numerical data. As it can be 
seen, the large distance behavior of pi is different in both 
cases, as for small values of IoidI/oo the OBDM presents 
an exponential decay, while in the opposite limit it is 
better reproduced by a power law. The comparison with 
the BHM shows qualitative agreement in the form of the 
decay, while quantitatively the description of the discrete 
model can be quite off, especially for strong interactions. 

The energy gap can be considered as an order parame¬ 
ter describing the insulating phase. Figure [^reports the 
gap calculated with different methods. In the first one, 
the (charged) gap is evaluated from the ground-state en¬ 
ergies calculated for A, A -I- 1 and A — 1 particles on 
L = A lattice sites, according to the expression 

Ac = Atv+i — 2E]\[ + Ejsf-i , (7) 

which corresponds to the difference of chemical poten¬ 
tials between the A -|- 1 and the A particle systems, re¬ 
spectively. Alternatively, an upper bound for the gap is 
obtained from the Feynman relation ([^. We obtain A 
and c by numerical fitting of the DMC data for S{k). The 
experimental data for the gap is taken from Ref. [M] . We 
observe a remarkable divergence between the exact DMC 
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FIG. 4: (Color online) OBDM pi{L/2) obtained from the 
pure DMC estimator for Vo = li?rec and two characteristic 
values of the interaction strength, |aiD|/ao = 0.1 (insnlating 
phase, left panel) and laiol/uo = 1 (snperfluid phase, right 
panel) for different system sizes L. Solid symbols, continuous 
model; open symbols, BHM. For |aiD|/(io = 0.1 one observes 
an exponential decay (insulating phase), while for |aiD|/oo = 
1 it is much better fit by a power law form (superfiuid phase). 



FIG. 5: (Color online) Energy gap. Squares with error bars, 
experimental data from Ref. [24] : filled circles, DMC predic¬ 
tion obtained from the fit to S{k —>■ 0); open circle, DMC 
prediction obtained according to Eq. Q; solid line, sine- 
Gordon model [23]; dashed line, upper bound given by Vo/2. 
uiD/uo = 0.18182. 


results, which are consistent among themselves using two 
different criteria, the SG model and the experiment. Still, 
the experimental points are in better agreement with the 
SG model than with the DMG calculation. This poses 
a question if the modulation spectroscopy method used 
in Ref. |21| is precise for measuring the value of the 
gap in shallow lattices. The gap Ac obtained from the 
BHM (not shown) grows from 0.9 for Vq — 0.5 to 1.5 
for Vo = 1.5 (in the units of i?rec) for the parameters of 
Fig. § and lies above the Vo/2 line. In this regime, Ac is 


larger than the energy gap between the Bloch bands and 
the BHM is not valid. 

To conclude, we established the zero-temperature 
phase diagram of a one-dimensional Bose gas in an op¬ 
tical lattice, determining the superfluid - Mott insulator 
transition line. We analyzed and compared the prop¬ 
erties of a continuous Hamiltonian (using DMC method) 
with that of the discrete Bose-Hubbard model (solved via 
exact diagonalization). We established the previously un¬ 
known regions of applicability of the approximate Bose- 
Hubbard and sine-Gordon models, and found that the 
sine-Gordon model fails to describe the regime of a shal¬ 
low lattice for any hnite lattice strength. This poses a 
natural question if it is possible reconcile the discrepancy 
by improving the sine-Gordon model. In general, Bose- 
Hubbard model is valid for high optical lattices with weak 
interactions, but the precise applicability of this descrip¬ 
tion depends on the quantity of interest. The depen¬ 
dence of the Luttinger parameter K on the height of the 
lattice and the strength of the interaction is reported. 
We also showed that the one-body correlations decays to 
zero following a slow power law in the superfluid phase, 
and exponentially in the Mott insulator phase. We com¬ 
pared our results with the experiment of Ref. |24] . and 
found an overall good agreement for the phase diagram. 
Instead, we saw a discrepancy in the value of the exci¬ 
tation gap. Importantly, our results help to understand 
the experiments with one-dimensional gases beyond the 
Bose-Hubbard approximation. 
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Note added .— After the present work was com¬ 
pleted and submitted for publication, a new experimental 
study |39] of the phase diagram by the LENS group in 
Florence appeared, in particular analyzing shallow lat¬ 
tices where we find discrepancy with the sine-Gordon 
model and the transport measurements of Ref. |24| . The 
experimental measurements and the path integral Monte 
Carlo calculations of Ref. [31] agree with our predictions. 
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